Operator normalized quantum arrival times in the presence of interactions 



G. C. Hegerfeldt and D. Seidel 

Institut fur Theoretische Physik, Universitat Gottingen, Tammannstr. 1, 37077 Gottingen, Germany 



Departamento de 



J. G. Muga and B. Navarro 
]utmica-Fisica, Universidad del Pats Vasco, Apartado Postal 644j 48080 Bilbao, Spain 



m 
O 
O 

> 
O 



> 
o 



m ■ 
o ■ 

^ : 

9 L,: 
^— > ■ 

G ■ 

3 : 



We model ideal arrival-time measurements for free quantum particles and for particles subject to 
an external interaction by means of a narrow and weak absorbing potential. This approach is related 
to the operational approach of measuring the first photon emitted from a two-level atom illuminated 
by a laser. By operator-normalizing the resulting time-of-arrival distribution, a distribution is 
obtained which for freely moving particles not only recovers the axiomatically derived distribution 
of Kijowski for states with purely positive momenta but is also applicable to general momentum 
components. For particles interacting with a square barrier the mean arrival time and corresponding 
"tunneling time" obtained at the transmission side of the barrier becomes independent of the barrier 
width (Hartman effect) for arbitrarily wide barriers, i.e., without the transition to the ultra-opaque, 
classical-like regime dominated by wave packet components above the barrier. 



PACS numbers: 03.65.Xp, 42.50.-p 



I. INTRODUCTION 

A major open question in quantum theory is how to 
include time observables and time measurements in a sat- 
isfactory way into the formalism. Pauli pointed out, in a 
famous footnote [1], that it is not possible to define self- 
adjoint time operators because of the semi-boundedness 
of the canonically conjugate Hamiltonian H (sec [2, 3] 
for recent discussions of Pauli's argument, its domain of 
applicability and implications) . A great deal of effort and 
ingenuity has been devoted to pinpoint and overcome this 
and other difficulties, in particular for describing the time 
of arrival of a quantum particle [3-39] . For an extensive 
review up to 2000 see [40]. 

In a physically motivated, axiomatic approach, Ki- 
jowski [4] derived an arrival-time distribution for free 
particles coming in from the left (or right), i.e. with only 
one sign of the momentum. For a one dimensional wave 
packet coming in from the left, with momentum (wave 
number) representation ip(k), k > 0, at time t = 0, and 
the arrival point at x — 0, Kijowski's distribution is given 
by 



njc(t) 



27rm 



dkijj(k)Vke- lhk2t/2m 



(1) 



Much more recently, this distribution has been related 
[12, 17, 23, 40] to the positive operator-valued measure 
(POVM) generated by the eigenstates of the Aharonov- 
Bohm (maximally symmetric) time-of-arrival operator 
[41], and generalized for systems with interaction [28, 31] 
and multi-particle states [33]. As pointed out in several 
reviews [2, 36, 40], one of the important pending ques- 
tions of Kijowski's distribution and its generalizations is 
the absence of an operational interpretation, in terms 
of some measuring procedure. This problem becomes 
particularly acute when interpreting the puzzling results 
obtained for complicated cases such as (i) freely moving 



particles with wave components coming from both sides 
to the arrival position, and (ii) particles subject to some 
interaction potential [34, 36]. 

Another research thread on quantum arrival times is 
based on modeling the irreversible detection by means of 
effective absorbing complex potentials. They were first 
used by Allcock in this context [5]. He chose an imag- 
inary potential step and found that a strong potential 
(therefore a short detection time) lead to important re- 
flection, i.e., many particles of the quantum ensemble 
associated with the incident wave packet could not be 
detected; this affected more intensely the slow end of 
the momentum distribution and distorted the measured 
signal. In the opposite limit, i.e., for a very weak ab- 
sorbing potential, the reflection was avoided but at the 
price of inducing a very large detection delay. By de- 
convolving the absorption-time probability density with 
the absorption-time distribution for a particle immersed 
at rest in the complex potential, in the limit of vanish- 
ing absorption, he could define an "ideal distribution" 
which turned out to be equal to the flux and therefore 
not positive definite, even for states formed only with 
positive momentum components. Kijowski's distribution 
for positive-momentum states emerged as the best posi- 
tively defined fit to this "ideal distribution" [5]. 

Twenty five years later, further investigations of the 
quantum arrival time using complex potentials revealed 
that perfect and fast detection were compatible by means 
of properly designed complex potential profiles [9, 22]. 
However, an ambiguity in the detailed form of these 
potentials remained, so it was not clear how to define, 
within the complex potential approach, an ideal time- 
of-arrival distribution in an arbitrary case. The distri- 
butions obtained with potentials constructed to absorb 
perfectly in the energy range of a given incident wave 
packet were in general close but not equal to Kijowski's 
distribution, and the relation with actual measurement 
procedures remained vague. A model proposed by Halli- 
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well of time-of-arrival detection by means of an unstable 
two-level detector made such a relation more plausible 
[19] (he derived a complex potential, again with a step 
form, from the original full set of equations for the two- 
level system), but the model lacked a precise physical 
content since the nature of the actual system or coupling 
mechanism were not provided. 

An important step forward to relate "ideal" and "op- 
erational" distributions has been the development of a 
model, using the "quantum jump" technique [42], for the 
detection of first photons emitted by atoms illuminated 
by a laser beam localized in a semi-infinite [35, 43] or fi- 
nite region of space [37]. Closely resembling Allcock's re- 
sults, it was found that strong driving by the laser caused 
reflection, whereas the weak driving regime could be used 
to get the flux by deconvolution from the first-photon dis- 
tribution, but Kijowski's distribution was not obtained in 
any limit [35]. The analogy with Allcock's work was not 
accidental since in the limit where the inverse life-time 7 
and the Rabi frequency Q become very large, but keeping 
the ratio fl 2 /j constant, and 7/fi >• 1, it is possible to 
find a closed equation for the atomic ground-state ampli- 
tude in terms of an effective complex potential [38, 43]. 
For the semi-infinite laser model, which is physically valid 
for slow enough atoms [37], and the laser tuned at res- 
onance with the atomic transition, it takes the simple 
form —ihilQ(x) 2 / (2-f), where Q(x) is the Heaviside func- 
tion, and the laser illuminates the region x > 0. In other 
words, once again, a purely imaginary absorbing step po- 
tential arises, but now with a clearcut physical content 
in terms of laser and atomic parameters. 

As pointed out before, though, the atom-laser model 
alone did not lead to Kijowski's distribution. In a very 
recent article, Kijowski's distribution has been obtained 
for freely moving states with positive momentum com- 
ponents [38], by applying the "operator normalization" 
proposed by Brunetti and Fredcnhagen [44] to the model 
that describes the interaction between the atom and a 
semi-infinite laser in the intuitive limit of strong laser 
field and fast decay (large SI and 7 respectively). 

From the previous discussion it should be clear that 
the problem of undetected atoms, i.e., atoms reflected 
(or transmitted in the finite beam case) without emitting 
any photon has to be faced in some way. For example in 
[35] "normalization by hand" for the fraction of detected 
atoms was used for strong driving conditions. Instead, 
for states with incident positive momenta the strategy 
of the "operator normalization" amounts to enhance the 
amplitude of each incident momentum component by an 
amount inversely proportional to the square root of the 
corresponding detection probability.Operator normaliza- 
tion may be applied whether or not the laser-atom pa- 
rameters allow for a simplified complex-potential descrip- 
tion, but the advantadge of such a parameter domain is 
that a much simpler one-channel calculation may be car- 
ried out; this is particularly useful to tackle the more 
complicated physical situations examined in the present 
work. 



One obvious limitation of the semi-infinite models 
(with the laser in explicit form or with a complex po- 
tential) is that it is not possible to study with them the 
arrival at a point, say x = 0, corresponding to a state 
incident from both sides, x > and x > 0, with nega- 
tive and positive momenta respectively; similarly, if one 
is interested in the arrivals in the midst of a potential 
interaction region, the semi-infinitely extended measure- 
ment will severely affect the dynamics of the unperturbed 
system on one side. The aim of this paper is to apply, in- 
stead of the semi-infinite interaction, a weak and narrow, 
minimally perturbing absorbing complex potential com- 
bined with "operator normalization" to these two more 
elusive arrival-time problems. 

The motivation for looking at the first case (free mo- 
tion) is a question originally posed by Leavens [16, 40]: 
When positive and negative momenta arc present in 
the wavefunction, Kijowski's original axiomatic approach 
does not apply and he obtained, in a heuristic way, a dis- 
tribution which in one dimension is given by a sum of 
independent contributions. The immediate consequence 
of this structure is that symmetrical and antisymmetri- 
cal states with respect to the arrival point x = 0, formed 
by adding or substracting two waves with opposite mo- 
menta, but identical otherwise, have the same arrival dis- 
tribution, i.e., the distribution is "blind" to the interfer- 
ence between the positive and negative momentum com- 
ponents and predicts arrivals for the antisymetric case 
even where the wave function vanishes at all times. This 
result may be understood formally from the resolution of 
the identity in terms of eigenstates of the time-of-arrival 
operator, and corresponds to a projection of the state 
onto the positive or negative momentum subspaces [40], 
but a hint on how that process could be implemented 
was missing [36]. We shall obtain this distribution by 
using the weak and narrow absorbing potential, i.e., an 
appropriate laser interaction, combined with "operator 
normalization" . 

The analysis of arrival times when the particle inter- 
acts with a potential has been problematic, too. It took 
a long time to generalize Kijowski's distribution in this 
case because the original axiomatic method is not appli- 
cable [13, 27, 28, 31, 45]. The proposed generalizations, 
though, are purely formal in character. In principle, a 
narrow and weak laser may be used to excite and detect 
the atom with minimal disturbance, even when the sys- 
tem's motion is affected by an additional interaction po- 
tential; by compensating the detection losses by "opera- 
tor normalization" a physically motivated generalization 
is obtained. As before, for simplicity, a complex poten- 
tial will play the role of the laser. In this paper we apply 
the method to tunneling across the paradigmatic square 
barrier potential, and obtain an arrival time distribution 
at the transmission side. It is similar to Kijowski's ex- 
pression, but modified by the phase of the transmission 
amplitude. The effect of operator normalization is that 
the mean arrival time becomes the average of the "phase 
times" with respect to the initial quantum state, instead 
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of the transmitted wave packet. As a consequence, the 
Hartman effect [46], namely, the essential independence 
of the average arrival time with respect to the barrier 
width, is obtained for arbitrarily large barriers, without 
the transition to an ultra-opaque, classical-like regime 
dominated by over-the-barricr momentum components 
[46, 48]. 



II. ABSORBING POTENTIALS AND 
QUANTUM ARRIVAL TIME DISTRIBUTION 

Let us consider an asymptotically free, moving wave 
packet impinging on an imaginary potential — iV e which 
is located in — e < x < e and a real potential U localized 
between a and b. The Hamiltonian is given by 



Hc = 2^, ~ iVtXe ^ + U X[a,b]( £ ) 



(2) 



where X[a,b] is one m [ a ib] and zero elsewhere, and 

/ _ / 1 : -e < x < e , , 



The scaling of the potential V e in the limit e 
by a function c(e) and 



V f = 



VqLq 
2c(e) ' 



is given 



(4) 



Here Vq and Lq are some arbitrary initial values for the 
potential height and potential width. Depending on the 
structure of c(e), we may investigate two cases: 

(a) c(e) = e, 

(b) c(e) ~ e a , < a < 1 . 

The case (a) yields in the limit e — > a delta-function 
potential with "strength" VoLq, V e Xe(x) — > Vo-£<o<K£)) 
whereas case (b) implies a weaker and less perturbing 
measurement. We shall eventually prefer this later case 
but include an analysis of the first one for completion and 
comparison. In the following we shall always apply the 

limit e — > and denote these two cases by e *-> 0, e 0, 
respectively. 

To obtain the time development under H c of a wave 
packet which is asymptotically free, one first solves the 
stationary equation 



H c 4>k = Ek<pk 
for scattering states with real energy 

h 2 k 2 



Ek 



2m 



(5) 



(6) 



By decomposing an initial state as a superposition of 
cigenfunctions of H c , its time development is obtained. 
This is easy for an initial free wave packet coming in from 
x = — oo in the remote past. Indeed, 



ip(x,t) 



pOC ____ 

Jo 



describes the time development of a state which in the re- 
mote past behaves like a free wave packet with ip(k), k > 
0, the momentum amplitude it would have at t = 0. (f>k 
would correspond in that case to the scattering states for 
left incidence. In the absence of a real potential (U = 0) 
we shall also be interested in states with both positive 
and negative momenta. The formal treatment for sym- 
metrical and antisymmetrical wave components is analo- 
gous, as shown later with more detail. 

It is convenient to use the interaction picture with re- 
spect to _ff — p 2 /2m, 



H I = e iH t/h( Hc _ Ho j e -iH t/h 
U 1 (t to) = e iH " t / h e ^ iH ^( t ^ t o)/ h e ^ iH " t »/ ri 



(8) 



where V\ is the conditional time development corre- 
sponding to H^. Then Eq. (7) can be written as 



^ = e-« ro */ ft C^(t,-oo)|^), 



(9) 



where \ip) = J dk\k)tp(k). An arrival time distribution in 
the region — e < x < e is given by the absorption rate 
[38] 



n(t) = -f J dx \^{x,tf 



(10) 



J dk dk' ^{k)^{k')e th{k2 - k ' 2)t/2m f e {k, k') 



(11) 



with the kernel function 



2V f e 

f e (k,k') = -^J dx <j> k {x)4> k i{x). (12) 

Now II(i) can be written as an expectation value on in- 
coming states in the form 



n(t) = (V|n t 



(13) 



and 



2V f 



fit - ~y u c(^ -™)hMu!(t, -oo). (w) 

To normalize on the level of operators we define [38, 54] 

/oo 
dt tl t = 1 - f/ c J (oo, -00)^/(00, -00) (15) 
-00 

and 

jjon = B-V2n t s-V2. (i 6 ) 

In Ref. [38] it has been shown that B is diagonal in k 
space, (k\B\k') — b(k,k')S(k — k') (remember that the 
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integrals over k run from to oo), and this leads to the 
normalized distribution 

n° N (t) = (mrw) 

= J dkdkjikjip(k r ) b{k,k)- L / 2 b(k' \k r )-^ 2 

ih(k 2 -k' 2 )t/2m 



x e" 



l fe(k,k'). 



(17) 



The kernel function b(k, k) can be calculated as in 
Ref. [38] or by the following simple argument. Because 
of / dt n ON 0) = 1, b(k,k)-^ 2 has to cancel the fac- 
tors which arise from integrating U(t) in Eq. (10) over t, 
namely 2 jjrf e (k, k). Therefore one has in the symmetric 
form provided by Eq. (16) 



b(k,k) 



-1/2 _ 



hk 



2-Kmf e (k, k) 



1/2 



(18) 





real potential 



U 



-iVe 



-£ +£ 



-£ +£ 



FIG. 1: Measurement scheme for the arrival time distribution 
at x — of a wave packet coming from the left, (a) Free 
wave packet, (b) Tunneling wave packet, (c) Arrival time 
distribution inside a potential. 



and 

n ON (t) 



2-7rm 



dk dk' ^{k)^{k')e lh{k2 - k ' 2)t/2m Vkk' 



fe(k,k') 



vOTa(P) 



(19) 



The normalization leads immediately to Kijowski's kernel 
\Jkk' in the arrival time distribution, modified by a model 
dependent term 



F e (k,k') = 



fe(k,k') 



y/Mk,k)f e (k',k') 



(20) 



which has to be investigated for specific situations in the 
limit e — > 0. 



III. EXAMPLES: FROM THE FREE CASE TO 
TUNNELING PARTICLES 

A. Arrival time distribution for free particles 

The easiest case to evaluate Eq. (19) is a free incom- 
ing wave packet of positive momenta. Its Hamiltonian is 
given by Eq. (2) with {7 = and the general solution of 
Eq. (5) in the region i, i = 0, 1, 2 associated with x < — e, 
—e< x < e, e < x, respectively, is given by 



i 



-iki'X 



(21) 



where fc = k 2 = k and k x = q £ = [k 2 + i^k^} 1 / 2 with 

lmq e > (Fig. la). The amplitudes Af and A~ are 
determined by the matching conditions at x — — e and 
x = e and they are explicitly derived in Appendix A. 

(i) Left incoming states: The appropriate eigenstates 
have boundary conditions = 1 and A 2 = 0. Then 



one can solve for the other amplitudes to obtain 



A+ = e~' ZlKt /D 



i ( k q t 



2 V* 



sm(2q e e)e- 2lk£ /D (23) 



-i(k+q e )t 



ID 



with the common denominator 



(24) 
(25) 



£) = cos(2g £ e)- + | ) sin (2(/,f). (2(i) 



In the limit e one has 



mV a L 



H 2 k + mVoLo 

h 2 k 



A + -> - 
2 h 2 k + mV a L 



Af, A 1 



h 2 k 



(27) 
(28) 
(29) 



2 h 2 k + mV L ' 

These results may be checked by considering a <5-potential 
in Eq. (2) from the start. 

(b) 



The limit e — > yields 



An 



o 

2 



A 

At,Al 



0, 

1, 
1 

2' 



(30) 
(31) 

(32) 



In this limit the measurement region has no effect on the 
motion of the atoms, as seen by Eq. (A10). However, by 
operator normalization a finite distribution is obtained in 
both limits. Inserting <jy£ \x) from Eq. (21) into Eq. (12) 
and calculating Eq. (20) leads to 



F e (k,k')^l, 6^0,6^0 



(33) 



■5 



and the arrival time distribution, (t), for a left inci- 
dent free wave packet equals Kijowski's distribution, 



iCT(t) - n K (t). 



(34) 



(it) The case of general states: The Hamiltonian H c in 
Eq. (2) commutes with the parity operator and, as a con- 
sequence, so does the operator Yl t in Eq. (14). Therefore 
its matrix elements between symmetric and antisymmet- 
ric states vanish so that the operator normalization of 
IT t can be performed in the subspaces of symmetric and 
antisymmetric states independently. 

Let us first consider an antisymmetric incident state 
\tp a ) composed of two identical wave packets with oppo- 
site momenta coming from the right and from the left 
[16, 40]. We first consider the antisymmetric subspace. 
The antisymmetric eigenstates \k;a) of Ho and 0£ of H c 
(k > 0) are given by 



\k;a) = -j=(\k)-\-k)) 



(35) 



and an antisymmetric wavefunction il> a can be decom- 
posed as 



IV>, 



a)= / 
JO 



dk |fc; a}(k; a\ip a ). 



(36) 



The antisymmetric eigenstates 4>t, of H c are obtained 
with the boundary conditions = 1 and = —1, 
which gives 

A+ = -Al 



A« = 



2 + *(|-f) sin ( 2( ^) 



-2ik< 



D 

— ike 



fcos(g e e) -ism(q e e)' 



(37) 
(38) 



The wave function in the presence of the imaginary po- 
tential can now be written as 



a (x, t)= dk (k; a\ij> a )<t>%{x) e ^ nkH l 2m . 
Jo 



(39) 



In spite of the fact that the wave function vanishes at 
x = the operator normalization preserves a finite arrival 
distribution even when the width of the measurement 
region contracts to 0. With Eq. (21) and Eq. (12) one 
has, using Eq. (20), 



F £ (M')-1, e^O.e^O, 



which yields 



27T77J 

h 

Tim 



dk (k;a\iP a )Vke- lhk2t / 2r ' 



/ dkii a {k)Vke- 1 
Jo 



Hk 2 t/2m 



(40) 

Vl) 
(42) 



In the last line we have changed to the ordinary mo- 
mentum representation and have taken the antisymmetry 
into account. 



A similar treatment may be applied to a symmetric 
wavefunction tp s by using symmetric eigenfunctions and 
the result is again of the form of Eq. (42) with ip a replaced 
by tp s - 

An arbitrary state 4>(k) can be written in terms of its 
symmetric and antisymmetric part as ip(k) — tp s (k) + 
tp a (k)- By parity, IL^ N (t) is the sum of the correspond- 
ing symmetric and antisymmetric contribution since the 
cross terms vanish. By means of a trivial calculation the 
sum can be written as 



ttON 



(*) = 



27T77J 



E 



dk^(±k)Vke- lhk2t / 2m 



(43) 

This is the operator-normalized arrival-time distribution 
for a general free wavefunction. The expression has been 
proposed in Refs. [4, 40] on more heuristic grounds as a 
generalization of the distribution in Eq. (1) from left (or 
right) incoming states to general free states. 



B. Arrival time distribution for tunneling particles 

A major advantage of the operational fluorescence 
model for the determination of arrival time distributions 
is that, in contrast to the approach of Kijowski, it is not 
restricted to free particles. This means that arbitrary 
potentials with bounded support can be considered and 
the arrival time distribution in the presence of these in- 
teractions can be calculated. 

For simplicity, we consider here the case of a rectangu- 
lar potential barrier. The Hamiltonian is given by Eq. (2) 
where, as before, the arrival at x — is measured and 
the additional real potential U is located in a < x < b. 
We here investigate the case with a < b < (Fig. lb). 

For solving the stationary Schrodinger equation, 
Eq. (5), the x-axis has to be divided into five regions 
i, i = 0, ...,4, corresponding to x < a, a < x < 
b,b < x < — e, — e < x < e, e < x, respectively. The 
general solution in region i is given by Eq. (21) with 
ko = k 2 = k 4 = k, ki = x = [k 2 - 2mU/h 2 } 1 / 2 and 
k 3 = q e = [k 2 + i^^L] 1 / 2 . In Appendix A we present 

the derivation of the Af using transfer matrices. 

In the case of an initial wave packet coming from the 
left and crossing the potential region the eigenstates re- 
quired have boundary conditions Aq = 1 and A J = 0. 
Then one can solve for A^ and A\ and obtain the ampli- 
tudes i.e., the solution in the measurement region. 

In the limit e one has with I = b — a and s = a + b 



ikl 



2 cos(>d) 

mVoLg 
h 2 k 



' I T + ^ ) ^ yJ ] 



+ e 



tks 



h k\ imVoLo 



(44) 
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whereas in the weak case the limit e 

p — ikl 

A ± - 

^3 



(*>) 



yields 



2cos(x0-*(f + ^)sin(xO' 



(45) 



which may also be obtained from Eq. (44) in the limit 
Vb — > 0. The last expression is independent of VqLq and 
of the position of the potential. It is half the transmission 
amplitude T(k) of a rectangular potential barrier without 
any arrival time measurement, see Fig. 2. Inserting Af 



R(k)e~ 



real potential 



U 



T(k)e 



ikx 
> 



FIG. 2: The stationary scattering function for a wave imping- 
ing from the left on a real potential barrier. T and R are the 
transmission and reflection amplitudes. 



into 



t>^ given by Eq. (21) and calculating F e {k,k') by 
Eq. (12) and Eq. (20) we obtain 



F e (k,k') = 



T{k)T{k') 
\T(k)\\T(k')\- 



e — > 0, 



(46) 



and the operator-normalized arrival time distribution of 
a tunneled particle takes, with Eq. (19), the form 



27T77J 



J dk i>{k)e- lRk2t/2m Vk 



T(k) 



\T(k)\ 



(47) 



The effect of the additional potential is the introduction 
of a phase factor, which is the phase of the transmission 

amplitude for the real potential in the weak limit e ^ 0. 

The dependence of the arrival time distribution H°g t (t) 
on the potential height U with potential width I = 10 is 
shown in Fig. 3. With decreasing U the arrival time at 
x = of a particle starting at (a;) = x < with mean 
velocity (v) — vq is first delayed but, for increasing poten- 
tial strength, an asymptotic distribution is reached with a 
mean arrival time ("Hartman time") tn = (\xq\ — 1)(v ) 
which is smaller than the free arrival time | xo | (v ~ 1 ) . This 
is related to the Hartman effect [46], which is discussed 
in detail in Section IV. The effect can indeed be observed 
[49] and it is complete in the limit U — ► oo, correspond- 
ing to an arrival point located in a "forbidden" region, in 
which case one has 



F e (k,k') -> e 



i(k-k')l 



U 



OO. 



(48) 



In the above setup an incoming free particle was pre- 
pared far to the left and then interacted with an ex- 
ternal real potential. A different situation arises when 
one includes the real potential as part of the prepa- 
ration procedure through which the particle passes far 



0.04 



0.03- 



^0.02 



0.01 



10 20 30 40 50 60 70 80 90 100 
t (a.u.) 

FIG. 3: Arrival time distribution II ON at x — in the pres- 
ence of a potential wall with height U = (solid), U = 0.3 
(dashed), U = 0.48 (dot-dashed), U = 0.58 (dot-dot-dashed), 
U = 1.0 (dotted), U = 2.0 (circles), potential width I = 10, for 
an initial Gaussian wave packet (xo = —50, Ax — 10, vo = 1) 
in atomic units m = H = 1. For U — > oo the mean arrival 
time approaches the Hartman time. 



away on the left and then continues to propagate freely. 
In this case the incident state for operator normaliza- 
tion purposes would be the normalized transmitted wave 
packet. Formally, this packet may be formed by a pro- 
jection onto a large region to the right of the poten- 
tial. For the positive results, i.e. transmissions, the 
normalized incoming free state is then characterized by 
T(k)i>(k)/(Jdk |T(fc)^(fc)| 2 ) 1 / 2 instead of ^(fc). Apply- 
ing Kijowski's distribution to the incoming free state thus 
prepared gives 

Il£(t) = J— I dki>(k)e" ihk2 ^ 2m VkT(k) 
zirm J 



J dk \T{k)^{k)[ 



(49) 



This expression coincides with the proposal of Refs. 
[13, 27, 31] for a generalized Kijowski distribution in the 
presence of an external potential and it is thus seen to 
be related to our result by a state preparation procedure 
that selects the transmitted particles. 

Remark: The transmission probability through a po- 
tential barrier as in Fig. 2 is given by / dk \^p{k)T(k)\ 2 
and one may argue that in this case the total arrival- 
time probability should equal this transmission probabil- 
ity. Instead of n°o t (£) one would then have a modified 

satisfying 



Instead of n°£ tV 
distribution, Tl°" t (t) 



J —oo J 



dk\^(k)T(k)\ 2 



(50) 



Injterms of operators this would require an operator, 
f[ t ON , satisfying 



r dttir= I 

J —oo J 



dk\T{kf 



(51) 
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With Eqs. (18), (20) and (46) it is easily seen that in this 
case the kernel of £? -1 / 2 takes the form 



70 r 



b(k,ky 

which yields 

F e (k,k')=f e (k,k') 



-1/2 _ 



hk 

2itm 1 



e — > 0. 



T(k)T(k'), 

The modified distribution n°Q t (t) is then given by 



(52) 



(53) 



ttON 
li pot 



(*) 



2nm 



J dk^(k)e- lhk2t/2m VkT(k) 2 , (54) 



which satisfies Eq. (50) and gives the joint probability 
density for both arrival and transmission. Normalizing 
this to 1 by hand just yields the distribution 11^ (i) of 
Eq. (49). Therefore 11% {t) can be understood as a con- 
ditional probability density for the arrival of the particle 
under the condition that it has been transmitted through 
the potential barrier. 




£45 



4QL 



0.5 1 1.5 

Barrier height U (a.u.) 



FIG. 4: The mean arrival time (t) of Eq. (55) behind a barrier 
of fixed width I — 10 (solid), I = 5 (dashed) as a function of 
the barrier height U. The initial wave packet is a minimal 
uncertainty Gaussian with xo = —50, vo = 1, Arc = 10 in 
atomic units m — h — 1. For C/ — > one has the free arrival 
time |a;o|(v _1 ), and for U — > oo the arrival time approaches 
the Hartman time tn = (|cco| — i)(u _1 ). 



IV. MEAN ARRIVAL TIMES 

To compare our result of tunneling times with previous 
works, it is useful to consider not only the arrival time 
distribution but also mean arrival times. We restrict our 
analysis to a wave packet that comes from the far left and 
collides with a rectangular potential barrier at a < x < b. 
The arrival time behind the barrier is measured at x = 
(Fig. lb). 

With T(k) = \T(k)\exp(i$ T (k)) the mean arrival 
time, (t) = J dttU ON (t), of the distribution in Eq. (47) 
is given by 



m f 



dk \i){k)\- 



[goj + & T (k) 
k 



(55) 



where xo < denotes the initial value for the mean po- 
sition of the wave packet. 

This result for (t) can be understood as the average of 
the "phase times" (the time required for a freely moving 
particle plus Wigner's time delay [50]) over the initial 
state. In contrast, previous proposals for the mean ar- 
rival time are written in terms of an average over the 
transmitted state [27, 46-48], which is just the first mo- 
ment of the arrival time distribution of Eq. (49). These 
results are not contradictory but correspond to different 
state preparations. 

The dependence with the potential height is shown in 
Fig. 4, where we have plotted (t) versus U for fixed barrier 
width I. In the free limit U — > 0, (t) approaches an 
"averaged free arrival time", |a;o|(w _1 ), and for U — > oo 
it approaches the "Hartman time" (\xo\ — ^(v -1 ) , where 
(v- 1 ) )= / dk^k^mk^/h. 

For analyzing the dependence with the barrier width / 
let us heuristically define a "mean tunneling time" as an 



easily calculable quantity obtained by subtracting from 
(t) the classical time for crossing the non-potential region 
with average momentum ho, 



to(|xq| - I) 
hko 



(56) 



The Hartman effect occurs when r in a tunneling colli- 
sion of a quantum particle with an opaque square barrier 
becomes essentially independent of the barrier width I 

This is shown in Fig. 5, were we have plotted the tun- 
neling time r versus the potential width I. For thin 
barriers r is above the free traversal time, as shown in 
the inset, but for increasing I there is a sudden transi- 
tion from a positive delay to a negative one for increas- 
ing U. It is clearly visible, that in the negative delay 
regime the tunneling time is nearly constant for increas- 
ing I (Hartman regime). In contrast we have plotted 
the tunneling time with respect to the arrival time dis- 
tribution n^(i) of Eq. (49), 7Tj which is related to the 
proposals of Rcf. [13, 27, 31]. Here the behavior of r for 
thin barriers is similar, but for increasing I the tunnel- 
ing time first slowly decreases and gets actually negative. 
This has caused many discussions and warnings that the 
interpretation of the "extrapolated phase times" as an 
actual tunneling time for the transmitted particle is un- 
justified [51, 52]. This question is not of our concern 
here though. The point we want to stress is that for very 
thick barriers the Hartman effect vanishes for tt, which 
for widths larger than a critical barrier length grows lin- 
early [47] . This is related to the influence of the exponen- 
tially decaying |T(fc)| in Eq. (49), which causes a dom- 
ination of the above-threshold components of the wave 
packet [27, 47]. By contrast, r, obtained with a different 
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preparation procedure does not show any transition to 
the classical-like, ultra-opaque regime. The intuitive ex- 
planation is that operator normalization compensates for 
all detection losses due to the l-dependence of \T(k)\, so 
that the result is never dominated by above-the-barrier 
components. 




Barrier width / (a.u.) 



FIG. 5: Mean tunneling time r (lines) of Eq. (55) versus I 
compared with the tunneling time tt (symbols) derived with 
Eq. (49) for different barrier heights U = (solid), U = 0.25 
(dashed, circles), U = 0.48 (dot-dashed, diamonds) and U = 
1.0 (dot-dot-dashed, triangles). The initial wave packet is a 
minimal uncertainty Gaussian with xo = —50, vo = 1, Ax = 
10 in atomic units m = H = 1. 



V. CONCLUSION 

We have shown that quantum arrival time distribu- 
tions can be defined in a physically motivated way by 
means of the absorption of a weak and narrow absorbing 
potential and operator normalization. As in Ref. [38] the 
use of the absorbing potential corresponds to the limit of 
an operational model, where a two-level atom impinges 
on a laser region and emits photons. By compensating re- 
flection and transmission losses in the absorption proba- 
bility by normalizing the corresponding operator function 
the arrival time distribution of Kijowski, IlR-(i), is recov- 
ered for freely moving states with positive momentum. 
The same method is applied to free states with arbitrary 
momentum components, and the result is a generaliza- 
tion of Hk (t) which agrees with the heuristic arguments 
of Ref. [4] and with the proposal of Ref. [40] . 

The weak and narrow absorbing potential allows a fur- 
ther generalization for interacting particles. For arrival in 
the transmission region behind a square barrier we have 
shown that our approach leads to an arrival time distribu- 
tion which only depends on the phase of the transmission 
coefficient. The mean arrival time is the average over the 
initial state of the "phase time" , so that the Hartman ef- 
fect, i.e., independence with respect to the barrier length, 
is obtained for arbitrarily wide barriers. 



A crucial ingredient for obtaining the present results 
is the way in which detection losses are compensated by 
operator normalization. Reasonable as this modification 
may be, there remains to be seen how this can be inter- 
preted from an operational point of view, i.e. how such a 
major intervention could be performed in practice. Oper- 
ator normalization can be viewed as a state modification, 
\ip) — > Z? _1 / 2 |^>), and we note that the modified state is 
not normalized. The physical process leading from the 
original state to the modified one may be a "filtering" 
preparation by scattering. However, the unbounded na- 
ture of the operator B~ x l 2 restricts such a method to 
states with a low-energy cut-off. The feasibility of the 
transformation in a more general case is an open ques- 
tion and further investigation is required. 
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APPENDIX A: TRANSFER MATRICES 

To calculate the transmission and reflection coefficients 
of the one dimensional potential problems of Fig. 1 or 
similar piecewise constant potentials we use the transfer 
matrix method of Ref. [53] and define "matching" matri- 
ces 

Mi(x) = ^ fc. e ik,x _f c . e -ik i x ^ • (Al) 

for each section i of constant potential Vi , where 

h = [k 2 -2mV l /h 2 ] 1 ' 2 1 (A2) 

with Imfci > 0, and fc, > for positive Vi. We shall al- 
ways use the index for the leftmost section, and number 
the rest rightwards and consecutively. We shall also de- 
note by Xi+i the boundary point between sections i and 
i + 1, so that the matching conditions take the general 
form 

M t (x t+1 ) = M, +1 (x i+1 ) . (A3) 

Multiplying both sides by M~ 1 (xi+i), 

where the transfer matrix connecting the amplitudes of 
regions i and i + 1 is given by 

T(i, i + 1) = Mr 1 (x i+1 )M i+1 (x i+1 ) (A5) 
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One may similarly obtain transfer matrices for non con- 
tiguous regions i and j > i + 2 by multiplying all the 
intermediate one-step transfer matrices, 

T(i,j) = T(i, i + l)T(t + + 2) T(j - (A6) 

In the following subsections we discuss the cases needed 
for the present paper. Note however that the transfer 
matrix method is of general applicability and may be 
adapted to continuous potentials leading in that case to 



differential equations for the amplitudes. 



1. Free particles 

For the free arrival time problem (Fig. la) one has 
three regions with fc = k 2 = k, k\ = q e , and match- 
ing points at x\ = —e and x 2 = e. From the matching 
matrices one obtains 



T(0,2) 



cos(2 9e e) - |(- 

i_( _k_ 
2\t 



2\q. 



+ f )sin(2g £ e) 
f )sin(2<z £ e) 



2ifee 



i(A_^ )sin( 2 gee ) 
cos(2g e e) + i(f + f)sin(2 ge e) 



(A7) 



and For waves incident from the left (x < 0) the bound- 

ary conditions arc Aq = 1 and A 2 = and there re- 
sults A+ = [T(0,2) n ]-\ Aq = T(0,2) 21 A+ A+ = 
T(l, 2) n A+ and A\ = T(l, 2) 21 A+ . 

(A8) 

This second matrix is useful to obtain the wave function 
in the absorbing region. In the limits e — I and e 
one has 




mvffio 1 h ^V„L B ) (A9) 2. Tunneling particles 

K*k 1 Wk / 

In the case of a tunneling particle (Fig. lb) with a < 

e ~*° b < one has to consider boundary conditions at X\ = 

lim T(l 2) = -( | (All) a ' X2 = b, %3 = — e and x 4 = e. T(2,4) is the same 

y 2 \1 1/ as the transfer matrix obtained in the previous section, 

^ /j j\ Eq. (A7), for the regions on both sides of the absorbing 

limT(l,2) = -f J. (A12) potential so we only need to calculate T(0, 2), 



T(0,2) = 



e lkl {cos(>cl) - |(| + ^)sin(xZ)} 



-|e" Jfcs (f - ^)sm(xl) 
- iH {cos(x/) + |(f + |)sin(xZ)} 



(A13) 



where I = b — a and s = a + 6. Since we are interested in 
the wave function inside the measurement region (— e, e), 
the corresponding amplitudes are given from the ampli- 



tudes of the rightmost region using the matrix in Eq. (A8) 
as before. 
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